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Abstract Colombeau's theory of generalised functions is used to calculate the contributions, at the 
rotation axis, to the distributional curvature for a time-dependent radiating cosmic string, and hence the 
mass per unit length of the string source. This mass per unit length is compared with the mass at null 
infinity, giving evidence for a global energy conservation law. 

o\ 
o\ 

1. Introduction 

>v 

An important problem in relativity is the calculation of distributional curvatures of non-regular metrics. 
Such curvatures represent physical fields whose support is confined to a proper submanifold of space-time, 
such as the point mass in the Schwarzschild solution and cosmic strings. 

The main obstacle to rigorously calculating such curvatures is the ambiguity that arises when one 
tries to multiply distributions, because the Riemann tensor is a nonlinear function of the metric and its 
derivatives. It was proved by Schwartz (1954) that one cannot define an associative multiplication of 
^ ■ distributions which is compatible with the usual operation of differentiation. There are metrics for which 

the distributional Riemann tensor may be calculated directly without the need to multiply distributions, 
but this class of metrics is very small. In fact a necessary condition for the metric to lie in this class is 
that the distributional curvature must have its support on a submanifold of codimcnsion of at most one 
(Geroch and Traschen, 1987); thus it is possible for metrics representing matter shells to lie in this class, 
f*** \ but not for strings and point particles. 

In an earlier paper (Clarke et al 1996) it was shown how one may apply Colombeau's theory of 
Generalised Functions (Colombeau, 1983, 1990), in which nonlinear operations are well defined, to the 
calculation of distributional curvatures where a direct calculations would not work. In Colombeau's 
5— i , theory, the Schwartzian space of distributions V(R n ) is embedded as a linear subspace of a much larger 

differential algebra G(R n ) by the smoothing operation 



- 1—^ 

X 



where $ G X>(R n ) has unit radius. Thus products of two distributions are defined in this algebra, 
although one cannot interpret this product as a distribution via the embedding. However one may give a 
distributional interpretation to this object via the notion of weak equivalence of two generalised functions 
F and G (written as F s=s G) which essentially is the generalisation of the equality of two distributions 
to G(R n ). 

In the previous paper (Clarke et al 1996), the focus was on calculating distributional curvatures of 
2-surfaces, in particular the cone with a metric of 

ds 2 = dr 2 +A 2 r 2 d(/) 2 , \A\<1. 

Here the Gaussian curvature density ^R^/g (as a generalised function) is weakly equivalent to the distri- 
bution 2n(l — A)S^(x, y), where (x, y) = (r cos</>, r, sin (/>). In this paper we shall consider applications to 
axisymmctric dynamical space-times with conical singularities and in particular space-times with a pure 
radiation energy-momentum tensor (Kramer et al 1980). That is with the energy-momentum tensor of 
the form 

T ab = 1> 2 k a k b , k a k a = 0, 

where vector k a is the null direction of propagation and $ represents the amplitude. In the case of 
cylindrical symmetry one may write the metric as 

ds 2 = e" 2 * (e 2 \-dt 2 +dr 2 )+ r 2 dcj> 2 ) + e 2 ^dz 2 , 



1 



where ip and 7 are functions of t and r. One may then obtain a pure radiation solution with 

kf, — kp — 1 

$ 2 = r~\ r - V/ - VV 2 
by solving the following field equations (Kramer et a_/ 1980; Krishna Rao, 1964) 

VVr + r"V,r - V'.tt = 
7,r + 7,* = r (^.r + "0,t) 2 

The Colombeau framework will be used to extract the distributional contributions at the axis r = 
for the Riemann tensor densities R al ' C d\f—Q and the corresponding distributional contributions for the 
energy-momentum tensor, which is defined via the field equations 

R a b — \Rg a b = T ab 

Such distributional contributions may then be used to give a definition of mass per unit length of such a 
string. 

Finally we shall consider the physical implications; the behaviour of matter at null infinity, using a 
three dimensional asymptotic formalism, developed be Ashtekar et al (1997). This will Demonstrate the 
sense in which global conservation holds. 

2. Distributional curvature of pure radiation fields 

In this paper we shall focus on the distributional curvature of a particular cylindrically symmetrical pure 
radiation solution (ip = 0, 7 = j(t — r)), where ^(u) is an arbitrary C°° function, whose metric may be 
expressed as 

ds 2 = e 2 7(t -r) ( _ d ;2 + dr 2 } + r 2 d(j) 2 + d ^ (1) 

We may write this metric in null Cartesian coordinates (u = t — r,x = r cos <j),y — r sin cj>, z); 

ds 2 = _ c 27(«) du 2 + 1 (dx 2 + dy2) + dz 2 _ 2c 2 7(M ) (l dudx + ^dudy) 

-\( X ^d^ + ^d X dy + y —^dyA. 
2 \ r z r z r z J 

Except at the singularity r = 0, the Riemann curvature density will have non-zero independent compo- 
nents 

R xy uxV^g= -i /r sin </>, R xy uyV^g = i/r cos <j>. 

which may of course be interpreted as distributions since such functions will always be locally integrable 
over R 2 . This does not however mean that this distributional curvature will be valid at the axis; we 
could expect a delta function contribution to the distributional curvature due to a non-zero angular 
deficit which may be shown to be 

lim27r|l- ^y*'" y /2 j = 2tt(1 - e-rW), 

as in the case of the static cone (See Clarke et al 1996). We shall be using a regularisation based 
upon Colombeau's generalised functions to recover the distributional contributions, on the axis, to the 
curvature. 

It may be remarked that this singularity is not quasi-regular but intermediate; that is the Riemann 
tensor components do not converge in the parallelly propagated orthonormal frame E along the integral 
curves of d/dr 

■.-id + .-*)! + *<i --">!: 



in which 



R 



0202 
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However there is an ortho-normal frame E' in which they do converge 



E' = 



e; 



d 1 — r 2 d 



+ 



1 + r 2 
2r L dt ' 2r 

I _ r 2 Q l+r 2 



2r 



dt 2r 



dr 



dr 



E'=^- 



d_ 

dz 



for which 



Rr 



0202 



- R' — R' — -re~ 27 V 
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We shall expect our distributional curvature to have cylindrical symmetry; this suggests that we 
should calculate g^ z by smoothing kernel &(u 7 x,y, z) e ,4i(R 4 ), with radius 



J \$(u,x,y, z)\ dudz > |, 



Ro=sup{(x 2 + y 2 ) 1 / 2 

The smoothed metric may then be written as 

&b e = -4 / 9ab(u' + u, x' + x, y' + y 1 z' + z)$(u'/e, x'/e, y'/e, z'/e) du' dx' dy' dz' . 

We shall expect our distributional curvature to have cylindrical symmetry; this suggests that it should 
be sufficient to use a cylindrically symmetrical smoothing kernel, which we shall denote as 

&(u, r) — <j>(u, r cos0, r sin0, z). 

With such a kernel the smoothed metric becomes 

cfob £ = ^ / 9ab{u' + u, x' + x, y' + y)®(u'/e, (x' 2 + y' 2 ) 1/2 /e) du' dx' dy' . 
£ Jr. 4 

The metric components will be in general sums of C°° terms, which may be identified with their 
smoothings in £m(R 4 ), and singular components. The only functions that we will have to smooth are 

(i) (u,x,y,z)^e 2i * 

(ii) (u,x,y,z) 1 — ► e 2 ~<( u )+^. 

For (i) the smoothing is given by (see Clarke et ai 1996) 



1 p 00 /-oc /-2-it j A(j> , / i0' \ 

? Li L e "*" ( "' /e - r ' /sy d -' d -' d *' 



ff,(r)e 2 '* 



where 



# e (r) = 4 



2tt 



00 /-r/e 



oo^O 



1 - 



e 2 r' 2 



$(u , r )r du dr , for r < e, 

2 /■OO r-Ro 

1 - 2tt— / / r' 3 <$>{u', r') du' dr', for r > e. 

r J-ooJO 



For (ii) the smoothing is given by 

e 2 7 («)+i0 = J_ / / e 2 7 («+«') 

6)3 J-ooJo JO 



;j.>\ V2 



re i<j>> _|_ r i e i4 



e i(0+0')/2 $ ( u '/ £j r ij e y du > dr > d(j) > 



in which (f)' may be integrated out by complex contour integration; If 

3 iW+*')/2 d 0', 



1 /2 

- / re i<t> + r > e i4>'\ 

o \ re 11 *' + r'e^ , 

then setting w = e 1 ^ and z = e 1 ^' , we may integrate out cj>' by integrating the complex function 

. / r'z + rw \ w 1 / 2 



r(z) = -i 

' rz + r'w J z 
around the circular contour k : \z\ = 1 to obtain 

f4B(r7r)e < * ) for r' < r, 
~ l4G(r/r')e i<?i , for r' > r 

where the hypergeometric functions E(\) and G(A) are defined by 



E(s] 



-an 



A 2 



dA 



2 \ V2 



A 2 S 2 

This enables us to express e^H^ as L e (u, r)e 27( " )+il ^, where 

/>r/e 



dA 



/oo /"rye 
/ c 2(7(MW) -T (u)) ^(er7r)$(u', r')r'du' dr' 
-ooJO 



L £ (u,r) = < 



OO /"OO 



3 2(7(u+£u')-7(«))/ 



ooJ r / e 
oo />oo 



Girjer , r )r du dr , for r < ei? , 



oo^O 



3 2( 7 (m+£m')~7(«)) ; 



E(er'/r)$(u', r')r' du' dr', for r > ei? - 



We may therefore express g a b e as 

oV = -c 2 ^ M W + I(dx 2 + dy 2 ) + dz 2 - 2e 2 ^L £ {u,r) ^dudx + ^dwdy) 



- 2 H s (r) 



9 9 9 9 

-i^-dx 2 + 4^da;dy + y _ dy 2 



The smoothed metric will be used to calculate distributional curvatures. In almost all cases, the following 
estimates for H e (r) and L E (u,r) will prove to be sufficient; 



Mr) = 
L e (u,r) = 



0(r 2 /e 2 ), for r < sRq, 

l + 0(e 2 /r 2 ), forr>ei? 

(9(r 2 /e 2 ), for r < sR , 

l + 0(e 2 /r 2 ), forr>ei? 



It must be emphasised that L e retains its dependence of u, even though we shall omit it from the order 
notation; on the other hand, H e is only dependent on r and e. 
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We now calculate the generalised Riemann density functions R ab c d V—g from the smoothed metric 



9ab £ 



R" 



R x \ 
R u 



-g e = Po + Pi cos 2(f> R uy uy ^g~ e = P - Pi cos 2</> 



R u 



-g e = R Uy u X eV Z Fe = PlS™2<l> 



R xy 



xv e 



-ge = -P'2 sin0 

^g~e = L E P 3 sin<f> ^ xy£ 

R Xy X y £ V^fe = P 3 



R uy , 



-ge = P2 COS (f) 
-g e = -L £ P 3 cos < 



where 



Po = 
Pi = 
P2 = 
Pi = 



0(r/e 2 ), for r < sRo, 
0(e 2 /r 3 ), for r > eR , 

0(r/e 2 ), for r < eRo, 
0(e 2 /r 3 ), for r > eR , 
—y(u)/r + 0(l/r), for r < eRo, 
-j(u)/r + 0{s 2 /r 3 ), for r > eRo, 

(0(l/e 2 ), for r <e.Ro, 
\0(e 2 /r 4 ), for r> eRo, 



To obtain the distributional parts to the curvature on the axis not contained in the expressions for 
R ab cd\/—g, one must evaluate the limits 



lim / ( R ab cd e - R ab cd V^9 j^dudx dy dz 



where \& is an arbitrary function in 2?(R 4 ). The following notations will be used throughout the analysis; 

K = supp 

K 1 = {{u,z) I {u,x,y,z)eK}, 

B e = { (u, x,y,z)€K\ (x 2 + y 2 ) 1 ' 2 < eR a } , 

R K = sup { (x 2 + y 2 fl 2 (u,x,y,z)e A"}. 

We shall first show that [R xy ux E y/^gl - R xy ux y/^g] w 0. That is 

J ^ ( R xy ux e V-fe - R xy uxV=g) * («, x, y , z) du dx dy dz -» 0. 



We have 



J ^ ( R XV ux e \T I ¥e -R XV ux * du dx dy dz 



< Mi 



Ib, 



R y ux e \l 9~e R y ux\/ ~H 



du dx dy dz 



-Mi 



L 



R xy ux-^f^g~R xy uxV~9 



du dx dy dz 



<M 2 e + M 3 -~ 
^ euo Kk 



where we shall use Mi to denote positive constants. Therefore 

[R xy uxeV^fe-R xy uxV^gj * {u, x, y , z) du dx dy dz = 0(e). 
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R 4 



Similarly it may be shown that 

[R ua u0 « 0. 



a, (3 = x, y 



For the component R ux xy \f—g (and similarly for R uy ' xy \/—g), we may write (using the mean value 
theorem) 

/ R ux x V £ v 7 -^* du dx d y dz = h + h 

JR 4 



where 



h = I R ux X y £ y/^f £ ^>(u,0,0,z)dudxdydz 

JK 

. 

h= I R ux 



-I 

JK 

[0,1] 



x Ve V -9e-^(u,£x, z)r du dx dy dz 



now 



Ib, 



R' J 



+ M 4 
(eR ) 



K—B e 
3 



xy £ V Se 
fyux 

n X y £ A 



r du dx dy dz 



r du dx dy dz 



< Mi 



1 1 



and 



/i = 



+ M 6 ^(— -— , 
eito Hk 



L £ P 3 rdr) / sin^d^ ) 0, 0, z) dudz 



IK' \Jo 

which vanishes because our radial integral has no ^-dependence and is bounded. This will imply that 



/ 

Jr. 4 



and hence that 



R v 



R 11 



xy 



-g £ ^ du dx dy dz = 0(e) 



-<j~0 a = x,y 



Only one component now remains; R xy xy \/—g, which we can expect to have a delta function 
contribution to the distributional curvature at the axis. We shall require a more accurate estimate for 
R xy xy e V~9e derived from 



R x \ 



-9e = -C 7 



dr 



This will imply that 

-Rk 



J 

Jo 



R xy 



x Ve 



-g £ r dr = e 7 1 



2tt f^ £ <f>(r')r' dr + l 

(i + H e y/*(i-H e + 2 ( ?~<L E 2 y 

2ir jf K,£ $(r')r' dr + l 



(1 + H e {R K )y/*{\ - H £ (R K ) + 2c27 Ls (u, R K f) 2 / 
But H £ (0) = L £ (u, 0) = 0, H £ {R K ) = 1 + 0(e 2 ) and L £ (u, R K ) = 1 + 0(e 2 ) so 



/■Rk 
Jo 



R xy 



xy e 



-g e rdr = e 7 (l-c- 7 ) + 0(e) 



We now consider the full four dimensional integral; 



/ 

Jr. 4 



R XV x V . V^T* du dx dy dz = h + I 4 



where for some £ € [0, 1], 



h 
h 



-L 



R xy xy E \/-^*(w,0,0, z) dudxdydz, 

R ^ xv e. \ Qe~ \ ( 



x Vs v 9s-^(t,^x,^y,z)rdudxdydz, 



Now 



\h\ < M 7 



R 11 xy e V 9e 



r du dx dy dz 



+ M 7 



< M> 



(sRq) 

8 2~ 



K-B, 
3 



R x \ 



'~9e 



+ M 9 e 2 f- 1 



r du dx dy dz 



1 



and 



Therefore, 



and so, 



\eR Rk , 
J (2^ {u) (1 - e~^ u) ) + 0(e)) *(u, 0, 0, z) du dz 

I R XV xy e yf^Fs^&udxdydz= [ 2TTC^ u \l~e-^)^(u 7 0,0 7 z)dudz + O{e) 

JR 4 JK> 

[R xy x V VH] ~ 27re^")(l - c-^)8^(x, y). 



The distributional contributions, at the axis, of the Riemann tensor may be used to calculate the 
distributional energy-momentum tensor density; 

[f^uy/Hl « T\^d- 27re^")(l - e^W)«( J )( a , y) 
» T Z Z V^ - 27re^")(l - e^)^*, y) 
P™6-\/-I] ~ T a bV^g (Other components) 



which still represents a pure radiation solution away from the axis, but whose distributional contributions 
on the axis are like those of a cosmic string. 

A physical interpretation to these contributions may be made in the sense of the mass per unit 
length of the string: Suppose that V a is a unit time- like vector in the direction of d/dt and that the 
distributional and non-distributional matter densities are defined as 

q = f ab V a V b 
g = T ab V a V b 

respectively. We may construct a generalised function of t and z, which we shall use to define the 
mass-per-unit length of the string via weak equivalence, by 



fi,(t, z) = 



d_ 

dz 



C(t,0,«) 



eV9 i3} - £V5 (3) 



dec' dy' dz' 



where C(t, a, b) is the constant t cylinder a < z < b of an arbitrary radius and g {3) and g (3> denote that 
the determinants of the induced metrics on the constant t hypersurfaces. 

Here V a = e<bl and 



l(3) 



-27 



1-H e + 2(2L £ - l)e 2 ^ 



This will imply that \/g (3) w e 2n '>/—g and hence, 



e~ 3 ^ u \g au f a u ^Tl~ gauT a u V^g 

_ e -7M ((ft. + £ £e f* u + y.L £ fy u )^qj 

r r J 



The first term will be weakly equivalent to 2tt{1 — c ' 1 ^)S^(x, y), where as the second term is weakly 
null, therefore on integrating we obtain 



H(u,z) « 2tt(1 - e- 7(u) ) 



(2) 



which implies that the string has a mass per unit length of 2tt{\ — c This is equal to the angular 

deficit since t = u at r = 0. 

3. Asymptopia and conservation of energy 

A natural question to ask is whether or not there is a conservation of energy law. At the local level one 
usually defines this concept in terms of the contracted Bianchi identities 



b;a 







which certainly hold if the space-time is at least C 3 -regular. If this was interpreted in the context of 
generalised functions, it would always hold since the smoothed metric is, by definition, C°°. 

In the case of radiating space-times, one may formulate the conservation law in a global sense; that 
is, whether or not there is mass loss between the source, in our case the string at r = 0, and null infinity 

In order to enable measurements to be made at null infinity, one must embed the space-time (M, g) 
into a compact space-time (M, g) whose boundary dM represents the points at infinity. In order to pre- 
serve causal structure under the embedding, one must impose the condition that g and g are conformally 
related; that is the existence of a smooth function on M such that 

g a b = Q 2 g a b 



and on dM, fl = but (3) V a O ^ 0. 

Many formalisms already exist for investigating asymptotic behaviour of four-dimensional space- 
times (Geroch, 1977). These are very often suited to axially-symmetric radiation solutions such as Bondi's 
radiating metric. A radiating string differs in that the singularity will itself extend out to infinity, thus 
making a choice of asymptote (M, g) problematic. In the case of our pure radiation solution, which is 
cylindrically symmetric, we are really measuring the mass-per-unit length of the string and can expect 
to measure the contribution of this piece of string at infinity. This suggests that a three-dimensional 
construction of null infinity would be most appropriate. 

Recently a formalism for constructing asymptotes at null infinity for cylindrically symmetrical radiat- 
ing space-times was devised by Ashtekar et al (1996) in which a symmetry reduction to a three-dimensional 
formalism is performed. This formalism was then used to investigate the asymptotic behaviour of Einstein- 
Rosen waves. Our application of this formalism to our pure radiation solution (ip = 0, 7 = j(t — r)) will 
be two- fold; firstly, to reveal to what extent it is asymptotically flat, and secondly to calculate the energy 
flux at null infinity so we can make comparisons with the matter measured on the string by distributional 
techniques. 
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3.1. Symmetry reduction 

The full four-dimensional metric (1) may be written in the form (Ashtekar et al 1996); 

ds 2 = e~ 2 ^ {3) g ab dx a dx b + e 2 ^dz 2 

where the 3-metric may be expressed in terms of Bondi coordinates (u = t — r, r, <j>) 

da 2 = (3) g ab dx a dx b = e 2 ~*(-du 2 - 2dudr) + r 2 d<t> 2 

The corresponding energy-momentum tensor for the 3-metric will be denoted by {3) T ab and may be 
calculated from Einstein tensor, via analogous three-dimensional field equations 

m Rab-\ m R m 9ab= {3> Tab 

In general {3) T ab will encode the metric coefficient ip as a scalar field. (See Ashtekar et al (1997) for 
Einstein- Rosen waves). In our case we have ip = so the only non-zero component of {3) T ab is 

(3) T _ 7(") 
J- uu 

r 

In order to examine the structure of null infinity (3) Z, we shall need to introduce a conformally related 
unphysical metric {3) g ab 

da 2 = fl 2 da 2 = e 27(M) {-f 2 du 2 + 2du df) + d(f> 2 (3) 

where Q, = f and f = r _1 . Here <3) X + is represented by the points f = 0. Since j(u) is arbitrary and 
smooth, we automatically have that {3) g a b is smooth across {3) I + . Moreover the hypersurface {3) I will be 
null; this is because we may define the normal vector h a = (3) g ab (3) Vf,fi, and 

^>g ab n a h b = f 2 e- 2 ^ ^0 (as r -» 0) 

3.2. Asymptotic flatness 

We now can define asymptotic flatness in the three-dimensional sense. The notation / = g will be used 
to denote that / and g are equal at the points of (3) X. 

Definition 1 (Ashtekar et al (1996)). A 3-dimensional space-time (M, (3) g) is said to be asymptotically 
flat at null infinity if there is an embedding into a manifold (M, (3, g) with a smooth metric and a boundary 
(3) X, topologically S'xR such that 

(i) M ~ {3) J is diffeomorphic to M. 

(ii) {3) g a b = Q 2 (3) g a b for some smooth function on M such that il = and (3) V a fi ^ 0. 

(iii) n (3) T ab is smooth and VL (3) T ab (3) g ab = 0. 

(iv) If h a = (3) g ab <3) V fc fi and V a is any smooth vector field on M, tangential to <3) I then fi" 1 i3) T ab n a V b = 


(v) If (3) V a n Q = then n a is complete on <3) X. 

Certainly conditions (i), (ii) and (v) are satisfied by our unphysical metric (3). Condition (iii) is also 
true since {3) T ab is trace-free and 

f! (3) T™ = -r 2 7» = 
However condition (iv) is not satisfied for our pure radiation solution since 

fT 1 {3) T ab h a V b = f- 1 (3) T UU ^g uf V u 
= i(u)e- 2 ^ u) V u 

The only way to make it vanish would be to force 7 to be a constant, thus making the string static. 
However the quantity f2 -1 {3) T ab n a V b will still be finite and smooth at {3) T. Condition (iv) was introduced 
in order to guarantee that the energy flux of matter across {3) X remains finite (Ashtekar et al 1996). We 
shall show that, despite this condition not being met, one may still define a mass at null infinity. 



9 



3.3. Bondi mass aspect and mass at null infinity 

Mass at <3> I may be defined from the energy-momentum tensors. Rather than working with the energy- 
momentum tensors {3) T ab and {3) T ab directly, we shall follow Geroch (1977) and Ashtekar et al (1996) and 
give definitions in terms of the 'modified' energy-momentum tensors 

(3 'L afc = n^s ab = nrn ab - \ <w> 5o6 ) 

(3> S ab = ^R ab -l {3) R l3) ga b 
The advantage of such tensors is that their algebraic relationship has the fewest number of terms; 

(3 >L ah = n (3) S ah + < 3) V a n b - IQ" 1 ^g ab ^g cd n c h d (4) 

We shall denote the pull backs to {3) X of (3> g a b, (3> L a b and {3) S ab as g ab , L ab and S ab respectively. 
The fact that the pull-back to {3) X of (3> \7 a n b vanishes and that (3) X is null will be sufficient to guarantee 
that g ab is a pull-back to (3) X of a positive definite metric on the space B of orbits of n a , which has S 1 
topology, and hence that g ab — m a m b , for some 1-form m a on <3) X. 

If it is the case that for any smooth vector field V a , i - 3) L ab n a V b = and (3) S ab n a V b = then Lah 
and Sgb have the form 

L ab = Lm a m b , Sah = Sm a m b 

There are sufficient conditions on the energy momentum tensor which enable {3) L ab h a V b = and 
{3) S ab h a V b = and are weaker than condition (iv) in definition 1. 

Proposition. Suppose that Til~ 2 , f = fl^ 2 {3) g ab n a n b are finite at {3) T and that for any V a , tangential 
to (3) J, {3) T ab n a V b = then the quantities (3) L ab h a V b and (3) S ab h a V b both vanish at (3 >X 

Proof. The vanishing of (3) L ab n a V b is clear from its definition 

(3) S ab = <3) T Qb -±T< 3 > 5ah 
=^ n- 1 i3) L ab h a V b = i3) T ab h a V b - \n- 2 Th a v b = (5) 
{3) L ab n a V b S 

To get {3) S ab n a V b to vanish, we should contract equation (4) with n a V b to get 

{3) Sa b n a V b = n- 1(3 'L ab n a V b - lfn a V a - i0t> o(3) V a / 
then apply (5) and the fact that / is finite at {3) I. 

With the pull-backs in the above forms, we arc able to define the Bondi mass aspect as 

Brh a = (S - L)rh a 

and if we denote a 1-dimensional cross section of (3) X, to which d/d<p is tangential, as C then the mass 
at (3) X for this cross section may be defined as 

E{C) = j>{\ - V2B~)m a dSa 

So with our solution we obtain a Bondi mass aspect of 

Bm a S* \e~^5i 

and the mass at <3> J for the cross section for u — constant (written as C u ) as 

E(C U ) = p (l - e-TM) deb 

Therefore the energy at infinity agrees with the mass per unit length at the axis (2), which implies that 
there is zero mass loss between the axis and null infinity. 
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